Detection of exchange interaction in STM measurements through Fano-like 
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We address Fano-like interference effects in scanning tunneling microscopy (STM) measurements 
of nanoscale systems, e.g. two-level systems. Common for these systems is that second order 
tunneling contributions give rise to interference effects that cause suppressed transmission through 
the system for certain energies. The suppressed transmission is measurable either in the differential 
conductance or in the bias voltage derivative thereof. 

PACS numbers: 73.23.-b, 72.25.-b, 68.37.Ef, 72.15.Qm 
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I. INTRODUCTION 

There are various techniques that allows one to de- 
tect and manipulate spin states in the solid state, which 
attract lots of interest. A partial list include opti- 
cal detection of electron spin resonance (ESR) in a 
single molecule^ tunneling through a quantum dot; 2 - 
and, more recently, ESR-scanning tunneling microscopy 
(ESR-STM) technique^ The interest in ESR-STM is 
due to the possibility of manipulating single spins^^ 7 - 
something which is crucial in spintronics and quantum 
information. Experimentally, modulation in the tun- 
neling current has been observed by STM using spin- 
unpolarized electron be&m££ Lately, there has also been 
a growing interest in using spin-polarized electron beam 
for direct detection of spin structures, 8 as well as uti- 
lizing the inelastic electron scanning tunneling spec- 
troscopy (IETS) for detection of local spatial variations 
in electron-boson coupling in molecular systems i^* 1 ^ 

Typically in STM measurements with an object lo- 
cated on a substrate surface, the tunneling current can 
either go directly between the STM tip and the substrate 
or go via the object. The tunneling electrons are thus 
branched between different pathways, which gives rise to 
interference effects when the partial waves merge into one 
in the tip or the substrate^ This interference leads to 
a suppressed transmission probability for the tunneling 
electrons at certain energies. The suppressed transmis- 
sion is a fingerprint of Fano resonances^ and generally 
appear in systems where tunneling electrons are branched 
between different pathways. Recently, Fano resonances 
have been studied in double and triple quantum dots 
systems ; 13 ' 14 ! 15 where the different pathways are consti- 
tuted of the different quantum dots. 

Fano resonances can be realized in a variety of sys- 
tem, ranging from systems with interactions between 
continuum states and a localized state, to systems where 
the branching of the wave function through diatomic 
molecules. The case originally discussed by Fano^ here 
reformulated for the purpose of STM and nanoscale sys- 
tems, interference occur between the different tunnel- 
ing paths in real space, one path going through the 



sample to the substrate whereas the other goes directly 
into the substrate. This idea is further exploited in 
works by Kubala and Konig,— Ladron de Guevara et 
aL»«, and Lu et al.J^- There are also many reports of 
IETS-STM measurements on molecules adsorbed onto 
metallic surfaces, which apart from the inelastic contri- 
butions also show Fanolike signatures in the transport 
dataj 16 ' 17 ' 18 ' 19 : 20 Theoretically, these features can be in- 
terpreted as interference effects between the different tun- 
neling channels . 2 1 1 2 2 1 23 1 24 

Here, we address the Fano interference effects that 
arise due to different pathways in phase space. In sin- 
gle level systems this phase space interference is can be 
caused by an effective spin-flip tunneling rate that is com- 
parable to the spin-preserving tunneling rate^ 5 - In two- 
level systems, tunneling paths such as | N — 2, n) — > \N — 
1, m) -> \N = 2, n') -»■ jJV = 1, m) -> \N = 2, n) give rise 
to phase space interference. Here, N = 1,2, denote the 
number of electrons in the state, whereas n, m are states 
indices^ Although such tunneling paths are of second 
order, they provide significant contributions to the trans- 
mission coefficient, and hence to the conductance, which 
accounts for the interference effect discussed here. The 
interference effects are sufficiently described in mean-field 
approximation of the sample correlation functions, hence, 
we do not discuss any fluctuations caused by electronic 
correlations or by the couplings to the tip and substrate. 

The paper is organized as follows. We introduce the 
formalism for the transport calculations in Sec. [Ill We 
discuss the Fano like interference effects and their impli- 
cations on two-level systems in Sees. IIII1 and we summa- 
rize the paper in Sec. IIVI 



II. TUNNELING CURRENT 

The STM system we have in mind can generally be 
described by the model Hamiltonian 



H = Htw 



H 



sub 



n, 



sample 



(i) 



where the first and second terms describe the electronic 
states in the tip and substrate, respectively. Here, we as- 
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sume flat band free electron like models for the states in 
the tip/substrate, and define H tip = J2 P aeti P £ p<? c l<T C p<7 
and H S ub = X^o-esub £ q^ c \a c q^i for the ti P and substrate, 
respectively. Henceforth, we let the wave vectors p (q) be- 
long to the tip (substrate). Creation (annihilation) of an 
electron at the energy e p{q)(J is denoted by c p{q)a (c p{q)a ), 
and we let a =f, \ denote the spin projection. The third 
term contains the information about the sample that is 
meant to be studied in the STM experiment, whereas 
Ht describes the tunneling between the sample and the 
tip and the substrate. The tunneling Hamiltonian can in 
general be written as 

Ht = ^ v pn<ycl a d na + ^ v qn<jc\ a d n<y + H.c, (2) 



where v p t q \ na is the tunneling rate between the tip (sub- 
strate) and the sample, whereas <v na (d na ) denotes cre- 
ation (annihilation) of a spin a electron at the nth level 
in the sample. 

The current is derived by standard methods, i.e. in the 
stationary regime we have J = J t i P = —ed/dt(N t i p (t)}, 
which give o 27 i 28 



J 



T(u)[fti p (u) - f suh (u)]duj, 



where the transmission coefficient 



T(u) = tr r tip G r Hr sub G a H. 



(3) 



(4) 



Here, G r ^ a - ) is the retarded (advanced) Green function 
(GF) of the sample which are to be determined later in 
the paper. The coupling between the sample and the tip 
is denoted by = 2%Y!,p v Zn<r v pmcr8(w - e p<7 ) while 

the coupling between the sample and the substrate is 
analogously defined by T^ a = 2TrJ2 q v* na v qma 5{uj - 
e qa ). In the expression for the current, Eq. ((3]), we 
have also defined the Fermi function /tip(sub) ( w ) = — 
/Utip(sub)) & t the tip chemical potential /iti p ( S ub)- F° r later 
reference we define the bias voltage eV = Hti p — ^sub, 
where we use utip = Mo + e V an d Msub = Mo with the 
equilibrium chemical potential /^o- 

The objective of the present paper is to describe the 
expected features in the transport characteristics that 
are caused by the Fanolike interference effects when the 
tunneling electron wave functions are branched between 
different pathways. It is thus sufficient to provide a mean- 
field expression for the GF of the sample, in which case 
the present level of transport equations, Eqs. and (TJ|, 
are valid. 



III. PROBING THE TWO-LEVEL SYSTEM 

Two level systems are the natural extension of sin- 
gle spins coupled to an environment of de-localized elec- 
trons. The most interesting physics in two-level systems 
is related to the lowest two-electron singlet state and the 



triplet states, and the exchange splitting between those 
states. Here, we will discuss an approach to measure the 
singlet-triplet exchange splitting parameter J in two-level 
system through the presence of Fano-like interference ef- 
fects. 

We begin our study by considering a diatomic molecule 
constituted by two identical atoms. One may think of 
the atoms as quantum dots (QDs) and the molecule as 
a double QD. We assume that the atoms are coupled 
through Coulomb and exchange interactions U and J, re- 
spectively. For simplicity, we assume infinite intra-level 
Coulomb interactions in order to avoid two electrons in 
one of the atoms, and we also assume that the tunneling 
between the atoms is negligible. This set of assumptions 
is not crucial for the effect we discuss, it merely permits a 
convenient framework for heuristic and qualitative stud- 
ies of the approach. The approach is straight forwardly 
generalized, which is done in the end of this section. 

The transport measurements are assumed to be per- 
formed by means of STM of a molecule located on a 
substrate surface. We therefore assume that the atoms 
couple equally strong to the surface, while the STM tip 
may couple asymmetrically to the atoms. We model the 
system by the Hamiltonian given in Eq. ([1]) where tip 
and substrate Hamiltonians are given in Sec. HH while 
the Hamiltonian for the molecule (sample) is given by 

"^sample = Zodna-dncr — 2 JSi • S2 

ray 

+{U- J/2)(m T + n u )(n 2T +n 2i ), (5) 

along with the condition that there may not be more 
than one electron in each level n = 1,2. Here, S„ = 
(1/2) Serer' ^L^'^u' ^ s a s P m operator, where a — 
(a x , a y ,a z ) is the vector of Pauli spin matrices. A com- 
plete model of the structure we consider may be given by 
Eq. (1) in Ref. [U 

In the given model, the sample can be either in the 
empty state |0), have one electron in either of the bonding 
or anti-bonding states 

IM(2)> = dim) ^ 2TU) |0), (6a) 

|i,3(4)> = < hm±^m \ Q)i (6 b) 

with corresponding energies E\ n = Eq, respectively, or 
have two electrons in the singlet 

. ^ ~ 4t d u ^ 

\/2 



2,1): 


= |5 = 


0,TO z = 


enerj 




= 2e + 


12,2) 


= \S 


= l,m z 


|2,3) 


= \s 


= z, m z 


|2,4) 


= \s 


= l,m z 



V2 



vt jt 



|0), (8a) 

(8b) 
(8c) 
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configurations, with energies i?2n = 2eo + U + J/2, n = 
2, 3, 4. Here we use the notation | N, n) where N = 0, 1, 2 
is the number of electrons, whereas n is a state label. 

The main effect that may be used for measurements 
of the singlet-triplet splitting J, arise due to phase space 
interference effects between the two-electron singlet and 
triplet states. This interference results in states that 
interacts only very weakly with the environment of de- 
localized electrons which, in turn, generate conductance 
suppression at biases that correspond to the singlet and 
triplet state energies. The conductance suppression is a 
direct response of the states that interact weakly with 
the surrounding electron bath. 

Qualitatively, one can understand the appearance of 
the localized states that interact weakly with the electron 
bath by considering the following. Assume that there are 
two electrons in the molecule and that they are config- 
ured in the singlet state |2, 1) = (dLdL — d\^d\,)\ff) / \/2. 
Removal of a spin f electron from the molecule can be 
done by removing an electron from either of the atoms, 
i.e. through the process (d^ + d 2 j)|2, 1) = — (d\, + 

4j.)|0)/\/2~ = -|1,4). Removal of a spin | electron 
thus always puts the sample into the anti-bonding one- 
electron spin J, state, which is orthogonal to the other 
one-electron states by construction. Analogously, remov- 
ing a spin I electron from the molecule being in the sin- 
glet state necessarily leads to a transition to the anti- 
bonding one-electron spin j state |1, 3). Clearly, the sys- 
tem cannot undergo first order transitions between bond- 
ing one-electron states |1, 1(2)) and the singlet state. In 
this respect one may therefore regard the singlet state as 
being decoupled from the bonding one-electron states. 

Likewise we find that the triplet states can be re- 
garded as being decoupled from the anti-bonding one- 
electron states, since removal of an electron from any of 
the triplet state configurations results in a bonding one- 
electron state. 

Having the sample described in terms of its eigenstates, 
it is beneficial to write the sample Hamiltonian in diag- 
onal form as 

^sample 

= J2 E Nn\N,n){N,n\. (9) 

Nn 

The tunneling Hamiltonian TLt is in the present con- 
text given by Eq. although we disregard spin-flip in 
the hybridization between the localized and de-localized 
electrons. Rewriting the operators d na in the tunneling 
Hamiltonian in terms of the eigenstates of the sample we 
obtain 

7~Lt = ^ ( ^ v pa N nm c pa + ^ v qa ]y nm c qa J 

Nnm * pa qa ' 

x\N,n)(N + l,m\+ H.c. (10) 
Here, the tunneling rates, c.f. Ref. HtJ 

V p (q)aNnm = {N, n\ {v p ( q ) la dla + V p ( q ) 2 ad2a) | N + 1, TTl) 

(11) 



also include the matrix elements for single electron tran- 
sitions in the sample. 

The sample GF in this system becomes a 20 x 20 ma- 
trix, which in general is not diagonal. This is because 
of higher order transitions that significantly contribute 
to the electronic structure and transport through the 
sample. The system can be simplified by assuming that 
only the diagonal processes like \N + l,m) — ► \N,n) — > 
| N + 1, m) and off-diagonal processes like \N + 1, m) — > 
| .AT, n) -> \N + 1, m!) -> \N, n) -> \N + 1, m) contribute 
to the tunneling. In fact, the off-diagonal processes we 
include into our scheme is sufficient in order to described 
the phase space interference in the present system. This 
simplification leads to a system of ten 2x2 matrix equa- 
tions, each of which is analytically solvable. 

Consider the GFs22£I 

i, in III (t,f) = H)(TA™ +1 (i)X™ + %(*')), (12) 

where AJ^ +1 = \N, n) (N + 1, m\. The discussed simpli- 
fications leads to the general equations of motion 

( i-Q^ ~ ^JVmn j Gjy nmm ' (t, t ) = S(t — t )P]SInmm' \p) 

Vn (t,t") 

-0 

X G Nnfi'm' 

(t",t')dt"(13) 

the transition energy AAr mn = E^+im — E^ ni the end- 
factors PNnmm'it) = ({^jvw+i' ^jv+WK*))i whereas 

VNnfj,fj,'n(t,t ) = ^""^ v paNnn v paNnn' 9pa(t, t ) 
pa 

+ V *qaNnpV q aNn t i'9qa{t, t'). (14) 

go- 
Here, gp( q )cr(t,t') is the GF for a free electron in 
the tip (substrate) satisfying the equation (id/dt — 

£p(q)a)9p(q)a(t, f) = S(t - t'). 

The GF in Eq. (fT3")) should be self-consistently solved 
for each value in the parameter space of the Hamiltonian 
Eq. {T]), bias voltage and temperature, according to pro- 
cedure lined out in Ref. HJ. Performing the calculations 
is such a way would provide a non-equilibrium description 
of the interference effects discussed here. For simplicity 
and in order to focus on the physical mechanism and 
effect we omit such a treatment. The self-consistent cal- 
culations are expected to change the quantitative results, 
however, the qualitative features will remain the same as 
in the present study. Here we discuss the physics with- 
out the self-consistency condition by further simplifying 
the equation of motion. Without losing information of 
the phase space interference effects, we assume that the 
end- factors PNnmm' — S mm ' ■ In this we assume that the 
off-diagonal occupation numbers are negligible, however, 
the off-diagonal GFs are non- vanishing and, furthermore, 
they provide important contributions to the widths of the 
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localized states in the sample. Moreover, we focus on the 
transitions between the one- and two-electron states. In 
absence of spin-flip processes the equations for transitions 
between the one-electron spin | and J, states are equal. 
It is hence sufficient to omit any reference to the spin 
degree of freedom, and therefore we consider processes 
that couple the two-electron singlet and states |2, 1) and 
|2, n), n = 2, 3, 4, through a one-electron state. 

It is worth pointing out that Eq. (fl~3|) is mean-held 
description of the sample, thus neglecting any type of 
fluctuations caused by the couplings between the sample 
and the electron baths in the tip and substrate. Fluc- 
tuations caused by the couplings to the tip and sub- 
strate electrons could by taken into to consideration by 
making use of the diagrammatic technique described in 
Refsj 30 i 31 . Especially interesting would be to also con- 
sider the singlet-triplet Kondo effect recently discussed 
by Paaske et al^ in the present context. However, a dis- 
cussion of the Kondo effect and other correlation effects 
is beyond the scope of the present paper. Here, we focus 
on the possibility of measuring the singlet-triplet split- 
ting J through direct transport measurements, which is 
sufficiently described within the present mean-field the- 
ory. 

Transitions between the different triplet state configu- 
rations do not give rise to the interference effects we dis- 
cuss in this paper. Therefore, it is sufficient to study the 
coupling between the singlet state and one of the triplet 
states at the time. It is, moreover, sufficient to con- 
sider only the interference between the singlet and, say, 
the triplet state |2,2). The interference arising between 
the singlet and the other triplet configurations merely 
renormalizes the coefficients in the final expression for 
the transmission. 

The Fourier transformed retarded GF for transitions 
between the two-electron states |2, 1) and |2,2) through 
the one-electron state n is then given by 



where 



UJ ■ 



A, 



Ai r 



jTjjl 



, (15) 



-2lmVl 
2^ 



ka 



C n (uj) 

mm'n (^0 
^1™^!""!'^^^^), (16) 



and T nm = T mnm , defines the combined coupling be- 
tween the sample and the tip and the substrate. Fur- 
ther, we define the coupling to the tip (substrate) by 



(q)a V p(q)alnm V p(q)<7ln™-' 



S(cu 



^tip 



Tmnm' • F° r a shorter nota- 



^tip(sub) 
1 mnm' 

such that r„ 

tion we have also put A rnn = Ai m „. Finally, the denom- 
inator C n (u>) — (lo — u> n +)(u> — contains the poles 
uj n ± of the GF. 



A. Symmetric coupling to the substrate 

Now, we assume that the atoms in the diatomic 
molecule couple equally strong to the substrate, i.e. 
Vqna = v qa . Then, the coupling matrices for the tun- 
neling between the sample and the substrate take the 
simple forms 



-psub 
1 |1,1> 



-psub 
1 1 1,3) 



sub 



r; 



-psub 
1 





1 

1 




(17a) 
(17b) 



We further assume that the tunneling rate between the 
tip and the molecule can be parametrized by 



J pna 



InVpa- The coupling matrices for the tunneling between 
the tip and the sample then become 



ptip 



■ptip 

11,3) 



r 



tip 



ntip 



(71 - I2) 2 -7? + 72 

-ll + 72 (7i + 72) 2 

(71+72) 2 -7i 2 +7 2 

-ll + ll (7i " 72) 2 



(18a) 
(18b) 



Finally, we assume that r* lp = Ar /2 and T s uh = T /2, 
where A <C 1, such that the broadening of the localized 
states T\2, m ) — ^ s \2 m>- This final assumption is made in 
order to simplify the analytical treatment and yields the 
GF poles loi± for transitions through the bonding state 



wi± =£ + Ut^- ~r (iTi), (19) 

and the poles u>3± for transition through the anti-bonding 
state 

w 3 ±=e + UT^-~T (l±l), (20) 

since A m i = A m 3, m = 1, ... ,4, and since the transi- 
tion energies An = E 2 \ — En = 2e + U — J/2 — e = 
eo + U-J/2 and A 2 i = E 22 -E n = 2e + U+J/2-e = 
Eq + U + J/2, which are associated with the singlet and 
triplet states, respectively. The poles of the GFs are 
thus located at the singlet and triplet energies, as ex- 
pected. Moreover, from these expressions it is legible 
that the poles lui + and cj 3 _ acquire vanishing widths, 
which marks energies for which the corresponding tran- 
sitions occur with vanishing probability. Hence, in the 
DOS of the sample there are states that interacts only 
weakly with the surrounding de-localized electrons. In 
this sense, there both broad and sharp peaks expected 
to appear at energies that are associated with both the 
singlet and triplet states. 

The sharply peaked state in the sample is a result of 
the interference between tunneling electron waves that 
are branched in phase space. This branching is respon- 
sible for the coupling of the singlet and triplet states 
through the one-electron states. The interference that 
arise between the branches of the tunneling electron 
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FIG. 1: (Color online) Typical DOS (bold) around the singlet 
and triplet states in the sample. The DOS associated with the 
poles u>i± (dashed, dotted) and a>3± (solid, dash-dotted) are 
plotted separately. Here we have the following parameters, 
£ = 2, U = 3, J = 0.6, To = 2 (units: meV), A = 0.005, 
7i = 0, and 72 = 1. 




4.2 4.6 5 5.4 5.8 

energy (meV) 

FIG. 2: (Color online) Transmission coefficient for the sys- 
tem with local DOS plotted in Fig. \T\ Total (bold) and 
partial transmissions through the channel which couples the 
two-electron state via the bonding (dotted) and anti-bonding 
(dashed) one-electron states. 



waves leads to both constructive and destructive inter- 
ference which manifest itself through the widths of the 
localized states in the sample. In terms of this argument, 
the states that correspond to the poles u>i + and 1^3- are 
subject to destructive interference which is sufficiently 
strong to remove any significant broadening of the state, 
it hence interact weakly with the surrounding electron 
bath(s). These states should display sharp peaks in the 
local DOS of the sample. The states that are associ- 
ated with the poles Wi_ and W3+ are subject to construc- 
tive interference which tend to increase the broadening of 
these states. Those states can, thus, be viewed as being 
strongly interacting with the environment. 

The plot in Fig. [T] (bold) illustrate a typical example 
of the calculated local DOS 

PtotH = ^2 Pnm (uj) = —^2,lmG\ nmm {u) (21) 



for energies around the singlet and triplet in the present 
context. There are two major broad densities which peak 
at the singlet (solid) and triplet (dotted) energies, respec- 
tively. These states, which have a strong interaction with 
the environment, mediate the tunneling current through 
the molecule. On top of the broad densities there are 
narrow peaks centered at the singlet (dashed) and triplet 
(dash-dotted) energies, respectively. These are the states 
that interact weakly with the de-localized electrons. In 
the calculations these poles have a finite width due to the 
small coupling to the electrons in the tip, c.f. the results 
for the poles given in Ref. [26|. The plots corresponds to 
the situation where the STM tip couples strongly to one 
of the atoms and only very weakly to the second. 

The states in the local DOS which interact weakly 
with the tunneling electrons are the source of dips in 
the differential conductance through the system. This 



can be directly seen in the transmission coefficient T{u>) 
for the system. Using the above assumptions and that 
the widths of the levels in the sample can be written 
as Tij = r*J p + 5ij^i2r| ub , we arrive at the transmission 
coefficient 



TM = 



(71 + 72! 
w - £ - U - J/2 



sq-U + J/2 



Ci{u) 



2 



(22) 



for the transmission through the two-electron states. 
This expression clearly shows that the transmission dips 
at the energies associated with the singlet and triplet 
states, and that the distance between the dips equals J, 
i.e. the exchange splitting parameter. 

The total transmission coefficient for the tunneling 
around the singlet and triplet states is shown in Fig. [2] 
(bold), along with the partial transmission coefficients 
for through channel that couples the two-electron states 
via the bonding (dotted) and anti-bonding (dashed) one- 
electron states. The plot demonstrate that the partial 
transmissions have an anti-resonance on the energy that 
corresponds to the singlet and triplet states, respectively. 
These anti-resonances are reflected in the total transmis- 
sion. 

Whenever the broadening of the quantum levels in the 
sample is larger than the singlet-triplet splitting J, there 
will not be two distinct peaks associated with the singlet 
and triplet state in the differential conductance. In addi- 
tion, visibility of the dips in the differential conductance 
require low temperatures since the transmission has to 
be convoluted with the thermal distribution functions, 
c.f. Eq. Therefore, in order to resolve the singlet- 

triplet splitting under those circumstances, one should 
preferably measure d 2 I/dV 2 rather than the differential 
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4.2 4.6 5 5.4 5.8 

bias voltage (mV) 

FIG. 3: (Color online) Transport characteristics for the sys- 
tem with local DOS given plotted in Fig. [T] a) Total dif- 
ferential conductance (bold) and the partial differential con- 
ductance through the channel which couples the two-electron 
states via the bonding (dotted) and anti-bonding (dashed) 
states, b) Total and partial bias voltage derivatives d 2 I/dV 2 . 
Here, T = 0.05 K, while other parameters as in Fig. [T] 



conductance. Those measurements provide further infor- 
mation about the long-lived states in the sample in terms 
of very sharp features at the energies corresponding to 
those states. 

In Fig. we plot the calculated transport character- 
istics corresponding to the DOS in Fig. [TJ The upper 
panel, Fig. [3] a), shows the total differential conductance 
dl/dV (bold) and the partial differential conductance 
for the tunneling through the channel which couples the 
two-electron states via the bonding (dotted) and anti- 
bonding (dashed) one-electron states, respectively. The 
broad electron densities give rise to wide conductance 
peaks at the singlet and triplet states which, however, 
are not clearly distinguishable in the total dl/dV since 
the level broadening is larger than the singlet-triplet ex- 
change splitting. Nonetheless, there are noticeable fea- 
tures, or dips, in the dl/dV at the energies for the sin- 
glet and triplet states. These dips are signatures of the 
states that are weakly interacting with the tunneling elec- 
trons. The large level broadening smears out the differ- 
ential conductance and makes a unambiguous read-out of 
the singlet and triplet states difficult. However, second 
derivative of the current d 2 I/dV 2 will provide the in- 
formation that is required for an unambiguous read-out 



since the conductance dips will display sharp features at 
the positions for the singlet and triplet states, as can be 
seen in Fig. [3] b). This panel shows the total (bold) 
and partial d 2 I/dV 2 for the tunneling through the chan- 
nel which couples the two-electron states via the bonding 
(dotted) and anti-bonding (dashed) one-electron states. 



B. Asymmetric coupling to the substrate 

Generally, the atoms in the molecule may couple asym- 
metrically to the substrate. The sample may, on the 
other hand, be comprised of a single entity with two lev- 
els. In both these cases we can assume that the hybridiza- 
tion can be parametrized according to v qn<7 = n n v q(T . 

The STM tip was already in the previous subsec- 
tion assumed to couple asymmetrically to the sample, 
hence, we merely have to re-derive the poles for the GFs 
and the transmission coefficients using the additional 
parametrization. The coupling matrices ^l^n) ^ a ^ e ^ e 
form of the ones given in Eq. (|18[) . however, with 7„ 
replaced by n n . We still assume that the tunneling be- 
tween the atoms in the molecule is negligible. Then, the 
localized states will be found at the energies 

wi± = s + U - ^T (kI + 4) (23a) 

4\/ (j- jKi« 2 r ) 2 - («; - «!) 2 (y) 2 

oj 3 ± = e + U- ^r («? + k%) (23b) 

4\Z( J+ 4 KlK2r °) 2 -( k "-^ 2 (y) 2 

for the transitions through the bonding and anti-bonding 
states, respectively. Again we assume that the level 
broadening is dominated by the coupling between the 
sample and the substrate surface. The expressions for 
the poles show that both solutions oj n ± , for both n = 1,3, 
acquire an appreciable width for asymmetric coupling be- 
tween the sample levels and the substrate. Thus, for suffi- 
ciently large asymmetry between the couplings, there will 
not be any states present that interacts weakly with the 
surrounding electron bath(s), and give rise to the conduc- 
tance dips. This is an indication that large asymmetric 
coupling between the levels in the sample and the sub- 
strate is detrimental to the ability of observing the sharp 
features in the dl/dV or d 2 I/dV 2 . Presence of such fea- 
tures will also depend on the asymmetry in the coupling 
between the sample levels and the STM tip. In general, 
in order to obtain localized states that interact weakly 
with the surrounding electron bath(s), the levels in the 
sample should be fairly symmetrically coupled either the 
tip or the substrate. Good symmetry in the coupling to 
both the tip and the sample enhances the amplitude of 
the sharp features in the d 2 I/dV 2 . 
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The transmission through the sample that is mediated 
by the two-electron states is modified according to 



T(ui) = A( ^ ) (71 «i + 7 2 K2^ 2 



u — £0 — U — qJ/2 



-E0-U + qJ/2 



C 3 («) 



(24) 



and in the sample become 



ptip 

A |M> 



-ptip 
± |1,3> 



ptip , 



1 1 



ptip . 

2 71 



-l 1 



(28a) 
(28b) 



where the asymmetry of the coupling introduces the fac- 
tor 



q 



7l K 2 + 72^1 

7i«i +I2K2' 



(25) 



Hence, apart from affecting the overall pre- factor, asym- 
metric coupling between of the levels in the sample and 
the substrate shifts the positions of the transmission dips. 
Quantitatively, coupling asymmetries 72/71 > 0.7 and 
K2/K1 > 0.7 yields about 6% shift of the transmission 
dips. Hence, even for asymmetric coupling this will be a 
reasonable good measurement of J. 



while the couplings matrices between electrons in the 
substrate and the sample are obtained by the replace- 
ments 7„ — * K n and tip — > sub. Here, we have only con- 
sidered the coupling matrices for the singlet state and the 
triplet state configuration |2,2) = \S = l,m z = 0). The 
coupling matrices that involve the other triplet configu- 
rations are the same. 

In the same approximation as previously, i.e. assuming 
that the broadening of the levels in the sample is domi- 
nated by the hybridization with the electrons in the sub- 
strate, we obtain the following poles to the GFs which in- 
volve the the two-electron states |2, 1) = \S = 0, m z = 0) 
and |2,2) 



C. Non- resonant levels 

The case thus far considered applies to system where 
the atoms are equivalent, so that the levels in the two 
atoms are resonant with one another. Generally, the 
two levels do not need to be resonant — they may also 
be spin-split, in which case the first term in the sample 
Hamiltonian is modified according to 



E 1 



no - d n(T d n(7 , 



(26) 



where e na — s n — erA n /2. Here, A„ is the spin-split of 
the nth level which is imposed by either internal and/or 
external magnetic fields. 

The singlet and triplet states remain the same as they 
are given Eqs. ([7]) and ([8]) while the energies for the 
respective state may change. The singlet state energy 
becomes E21 = £1 + £2 + U — J/2, whereas the triplet 
state energies become E22 = £i+£ 2 + £7 + J/2, E23 = 
E22 - (Ai + A 2 )/2, and E 24 = E 22 + (Ai + A 2 )/2. The 
one-electron states is in this modified model preferably 
described by 

\l,n) = {<5„i<4 t +<*n2<4i +5„s4t + 6 n4dl l }\0), (27) 

with corresponding energies E ln = (S n i+S n3 )e n ^ + (S n2 + 

We consider the general case in the sense that we relax 
the assumption of the levels being symmetrically coupled 
to the substrate. We thus assume that the tunneling 
rates can be parametrized according to v pna = ^nVpa 
and Vq na = K n v qa , as in the previous subsection. The 
coupling matrices between the electrons in the STM tip 



UJi± 



^3± 



£i| + £2| + 2U 



^2 r 



±2 V J2 - («2lV2) 2 

= - 2 7«ir 

1 



±^/J 2 -(nlr /2y 



(29a) 



(29b) 



The poles u>2± and lo^± , for the spin J, channels, are ob- 
tained by letting ]— >{ in the equations above. From the 
form of the poles we deduce that, in the spin-degenerate 
case there may be sharp localized state only in the case 
when J -C K„ro/2- In the symmetric case it should be 
possible to measure J in cases where the singlet-triplet 
splitting is small. Then, the theory presented in the pre- 
vious subsections apply, and therefore we proceed our 
discussion for the spin-dependent case. 

The spin-dependent case we consider is based on the 
assumption £1 = £2 = £0 and that the level spin-split 
is uniform, i.e. Ai = A 2 = Ao. We also assume that 
the levels in the sample couple symmetrically to the sub- 
strate. In this case we can use one-electron states given 
in Eq. ^ and the coupling matrices given in Eqs. (|17p 
and (|18|) . Considering spin-dependent transport enable 
studies of the spin-splitting of the two-electron triplet 
states. The singlet state couples to all triplet configura- 
tions, through the one-electron state, which therefore will 
display sharp localized state at different energies. These 
different energies correspond to the spin-splitting of the 
triplet configurations. Due to the spin-dependence, there 
appear four sets of poles which couple the singlet with 
each of the triplet state through the one-electron states. 
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For transitions to the one-electron spin f state we have 



IV. SUMMARY 



2,2) 

A± 

,|2,2) _ 



^ + ^ + C/T^-^r (l T l), (30a) 



V = £o + ^ + C/T^-^r (l±l), (30b) 



J i. 



2 8 



12,3) 12,2) a / -« -i \ 

= - A (l T 1), 



.12,3) 
^3± 



|2,2 



r -Ao(lTl), 



(30c) 
(30d) 



while the poles associated with transitions to the one- 
electron spin J. states are given by 



,|2,2) 



A 



J i. 



± = eo-^ + UT g - gro(lTl), (31a) 
4± 2> = £ o-^ + C/T^-^ro(l±l), (31b) 



.12,4) 
J 2± 



|2,2) 



+A (lTl), 



12.4) 12.2) a / -. 1N 



(31c) 
(31d) 



Note that the poles are equal in pairs according to 

Due to these 



12,2(3)) 

J i± 



12,4(2)) , |2,2(3)) 



^3± 



12,4(2)) 
J 4± 



equalities we expect four sharp peaks in the local DOS 
of the sample which are associated with the two-electron 
states, two for the singlet and two for the triplet states. 
The reason is clear since the one-electron spin f states 
do not at all couple to the triplet state |2, 4) = dLdJ, |0). 
Therefore, there cannot appear any sharp peak associ- 
ated with the corresponding transitions. The same argu- 
ment applies to the spin j channel. In any other respect, 
each spin channel can be treated separately by means of 
the theory developed in Sees. IIII Al and IIIIBI 



We have presented a theoretical prediction of the pos- 
sibility to measure the singlet-triplet exchange interac- 
tion parameter J through Fano-like interference effects. 
We argue for measurements of d 2 I /dV 2 for extraction 
of J, since the level broadening in many realistic situa- 
tions is larger than the singlet-triplet splitting. The level 
broadening therefore effectively smears out possible iden- 
tification features in the differential conductance which 
makes the read-out of J difficult. We also address the is- 
sue about asymmetric coupling between the sample levels 
and the substrate, and find that moderate asymmetries 
preserve a reasonable good measurability of J. Finally, 
we address the question of non-resonant levels. There will 
not be any Fano-like interference effects between which 
enable singlet-triplet splitting read-out unless the levels 
are resonant. However, uniform spin-splitting of the lev- 
els do not destroy the measurability of the singlet-triplet 
read-out. In the spin-split case we predict four sharp fea- 
tures in the d 2 I/dV 2 instead of two, which enable pair- 
wise read-out of i) the spin-split and ii) the singlet-triplet 
splitting. Experimental results on two-level systems us- 
ing the addressed set-up would be very intriguing and 
may open novel approaches to information storage. 
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